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We investigate neutrino oscillations in constant matter within the context of the standard three 
neutrino scenario. We derive an exact and simple formula for the oscillation probability applicable 
to all channels. In the standard parametrization, the probability for v e — » transition can be 



written in the form P(y £ 



A S fj, cos 5 + B e fj, sin S + C eM without any approximation using CP 



phase S. For — > v T transition, the linear term of cos 28 is added and the probability can be written 
in the form P(y^ — » u T ) = A^ T cos 8 + _B MT sin 8 + (7 pT + D^ LT cos 28. We give the CP dependences of 
the probability for other channels. We show that the probability for each channel in matter has the 
same form with respect to 5 as in vacuum. It means that matter effects just modify the coefficients 
A, B, C and D. We also give the exact expression of the coefficients for each channel. Furthermore, 
we show that our results with respect to CP dependences are reproduced from the effective mixing 
angles and the effective CP phase calculated by Zaglauer and Schwarzer. Through the calculation, 
a new identity is obtained by dividing the Naumov-Harrison-Scott identity by the Toshev identity. 

PACS numbers: 14.60.Pq 
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I. INTRODUCTION 

Strong evidences for neutrino oscillations have been ac- 
cumulated and they indicate the finite masses and mix- 
ings of neutrinos. The first evidence is the observation 
of zenith-angle dependence of the atmospheric neutrino 
deficit [[l], which is consistent with — > v T transition 
with the mass difference and the mixing as 



A. 



32 



3 x l(T 3 eV 2 , 



sin 2 29 



23 



1. 



(1) 



where = raf — mj and rrii is neutrino mass. The 
second evidence is the solar neutrino deficit 0, which is 
consistent with v e — > v^/vt transition. Recent reports for 
the measurements of CC and NC interactions in the SNO 
experiment || are consistent with the Standard Solar 
Model and strongly suggest the LMA MSW solution 



A 



21 



7 x 10 _b eV^ 



sin 2 20 



12 



o.. 



(2) 



to the solar neutrino problem. Thus, the information 
on neutrino masses and mixings (Maki-Nakagawa-Sakata 
matrix) [Q has been increased considerably. The main 
physics goals in future experiments are the determina- 
tion of the unknown parameters #13 (only upper bound 
sin 2 20i3 < 0.1 is obtained S) and the CP phase 8 in 
addition to the high precision measurement of the pa- 
rameters in (ffl) and (||). In particular, the observation 
of 6 is quite interesting from the point of view that 6 is 
related to the origin of the matter in the universe. In 
order to get at the CP phase 5, a considerable number of 
studies have been made on CP violation H p|, pj, p|, fL0|. 
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Moreover several experiments using neutrino beam are 
planned JTT| , p"2| . In such experiments with high energy 
neutrino beam, matter effects cannot be neglected and 
prevent the precise measurement of pure CP effects. In 
this paper we derive an exact and simple formula of os- 
cillation probabilities for all channels. 

Before giving our results, let us give the present sta- 
tus of exact formulas associated with oscillation prob- 
ability in constant matter. The mass eigenvalues have 
been given by Barger et al. [ p"3] . Zaglauer and Schwarzer 
]l4| have calculated effective mixing angles and effective 
CP phase in the standard parametrization. Furthermore, 
the elements of the MNS matrix have been calculated 
by Xing fl5[| in parametrization independent way. How- 
ever, these exact formulas obtained by diagonalizing the 
Hamiltonian, are complicated and it is not easy to ob- 
tain useful information on matter effects and the CP 
phase. Ohlsson and Snellman have directly calcu- 
lated the amplitude without diagonalizing the Hamilto- 
nian although the CP phase is not considered (which is 
similar to our result as shown in Sec. III). 

Recently, by comparing the Hamiltonian in matter and 
that in vacuum, Naumov |l7j and Harrison and Scott |lg}| 
have presented simple and elegant identity 



A12A23A31 J = A12A23A31 J, 



(3) 



where J is the Jarlskog factor |l9| and tilde represents 
the quantities in matter. We have investigated the mat- 
ter enhancement of J in Ref. ^(J by using this identity. 
Furthermore, the enhancement of T- violating quantity 
APt = P{v e — * Vfj) — P{Vfj, — * v e ) has been also inves- 
tigated by Parke and Weiler [^J and Akhmedov et al. 

In our previous paper |p3) , we have found other identi- 
ties on the product of the MNS matrix elements by com- 
paring the Hamiltonian in matter and that in vacuum. 
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From these identities, we have presented the exact for- 
mula for P(u e — * Vn) and we have definitely separated 
pure CP effects from matter effects. In the standard 
parametrization, we have also shown that the probability 
for v P — > v n transition can be written in the form |331 P3l 



A eil cos 8 



B efl sin S 



C e 



(4) 



It has also been found that CP trajectory proposed by 
Minakata and Nunokawa |^4j is exactly elliptic in bi- 
probability space when 5 changes from to 27r. 

In this paper, extending our previous results ]2!j[ ], we 
derive an exact and simple formula of the oscillation 
probabilities for all channels. We also find that the prob- 
ability for each channel in matter has the same form with 
respect to the CP phase 8 as in vacuum. For example, 



the probability for 



transition can be written as 



P{ v n - * v t ) = A^ T cos 5 + B^r sin 5 + C^r + D^ T cos 26, 

both in vacuum and in matter. Namely, matter effects 
just modify A^ T , B^ T , C^ T and D^ T . We also give the 
exact expression of these coefficients. Furthermore, we 
show that our results with respect to CP dependences are 
reproduced from the effective mixing angles and effective 
CP phase calculated by Zaglauer and Schwarzer jljj . Al- 
though the effective CP phase are complicated, the CP 
dependence of probability is simplified for all channels. 
This is guaranteed by the Toshev identity 

S23C23 sin 6 = S23C23 sin 8. (6) 

Finally, we obtain a new identity 

A12A23A31S12C12S13C13 = A12A23A31S12C12S13C13, (7) 

related to the 1-2 mixing and the 1-3 mixing by divid- 
ing the Naumov-Harrison-Scott identity by the Toshev 
identity. 



II. CP DEPENDENCES OF OSCILLATION 
PROBABILITIES IN VACUUM 

In this section, we review the CP dependences of the 
oscillation probabilities in vacuum. The Hamiltonian in 



H = 



is diagonalized as 



Hp p Hp .11 Hp 



H p e H^^ 



Hftr 

Hre H T n H TT 



(8) 



U*HU = 




A 2 i 



A 3 i 



(9) 



where U is the MNS matrix which connects flavor eigen- 
states with mass eigenstates [pi] . Then, the amplitude 
of v a to vp transition after transporting the distance L 
is 



A{v a -> up) = J2 U< 



* e-t^Upi 



(10) 



i=l 



where a and (3 denote flavor indices e, /1 or t. And the 
oscillation probability is given by 

P[y a -> up) = \A(u a -> u p )\ 2 

cyclic cyclic 

= S af3 - 4 £ ReJ^ sin 2 A^ ± 2 £ Jsin2A^.,(ll) 



where 



J^ = U ai U^{U aj U^)\ J = lmJ^, (12) 



AL 



A,jL _ (to 2 - m 2 )L 



4E 



4E 



(13) 



The ± sign of the third term takes — (+) in the case that 
(a, /3) is given by (anti) cyclic permutation of (e, /i). The 
cyclic sum is over (ij) = (12), (23), (31). 

Next let us calculate Re Pj^ and J to study the CP de- 
pendences of probabilities in the standard parametriza- 
tion S, 



where Sy = sin dij , Cy 
Re Pj^ are given by 



-S12C23 
S12S23 - 



C12C13 

Cl2C23Sl3e lS 



ReJ, 
ReJ, 
ReJ. 



12 

23 

efj, 

31 



C12C23 " 
-C12S23 



S12C13 
- Sl2S23Sl3e l5 
Sl2C23Sl3e Z 




-(c 



J r COS 8 — S19S 



\ 2 )Jr COS 8 + 
2 „2 2 2 



2 2 



2/2 2 _ 2 \ 
12 c 12 13l s 23 s 13 c 23/ 



I2 a 23 a 13'-13) 

t r 2 2 2 2 

-J r COS — C 12 S23 s 13 c 13i 



(14) 



(15) 
(16) 
(17) 
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and ReJll are given by 



ReJ^ = {c 2 2 - S 2 2 )J r COS 8 + S 2 2 C 2 2 C 2 3 (c 2 3 S 2 3 -S 2 23 ), (18) 

ReJ 2 ? = - J r cos 8 - s 2 2 c 2 23 s 2 13 c 2 13 , (19) 
ReJ^ = J r cos (5 - c\ 2 c 23 s\ 3 c\ 3 , (20) 



and ReJ z J T are given by 



r> t12 2 2 2 2 /-|,2,4\_/2 2,2 2 \ 2 

" eJ /iT — S 12 c 12 s 23 c 23l i ' s 13 "+~ S 13/ V s 12 c 12 + s 23 c 237 s 13 



( C 12 - S 12)( C 23 _ S23) S 12 c 12«23C23Sl3(l + «i 3 ) COS S + 2s\ 2 <? 12 S 2 23 (? 23 s\ 3 COS 28, (21) 
'/.;- ~( c 23 ~~ S 23) Jr COS 5 + S 23 C 23 C 13 (s 12 S 13 — C 12 ), (22) 

ReJ^ = (c2 3 -S2 3 ) J rCOS(5+S23 C 23 C 13( C 12Sl3- S 12), (23) 



where J r = S12C12S23C23S13C13. We also obtain the Jarl- 
skog factor, 



J = J r sin 8. 



(24) 



It should be noted that the replacements S23 — * c 2 3, 
C23 — * -S23 in ReJ*^ leads to ReJ^. 

In summary, the CP dependences of transition prob- 
abilities are obtained by substituting ([l5j)-(|4]) into (|T~ 

as 



A eil cos 5 + B sin 8 + C eM , 



P(v e — > v T ) = —A eli cos 8 — B sin 5 + C eT , 



(25) 
(26) 



cos 8 + B sin S + C^ T + D cos 25.(27) 



P{v e — > is ft) and P(y e — > v T ) have only linear terms of 
cos 8 and sin 8. On the other hand, Piv^ — > v T ) has also 
a term proportional to cos 28 in addition to the terms of 
cos 8 and sin 8. 

Next, let us consider survival probabilities. One can 
calculate the CP dependences from the unitarity and the 
transition probabilities. Those are given by 



P{Ve 

P(v» 



C e 



(28) 



cos 8 + C w - D cos 28, (29) 



P{v T — > v T ) = A TT cos 8 -\- Ctt — P) cos 28, (30) 



where the coefficients can be calculated from (|l5|)-(24). 
Note that the term proportional to sin<5 disappears in 
survival probabilities. In particular, not only sin<5 term 
but also cos 8 term does not exist in P{v e — > v e ) |28|]. In 
P(Vfn — > Vfj) and P(y T — > v r ), the linear term of cos 28 
has opposite sign to that in P{v^ — » i/ r ). The proba- 
bilities for CP conjugate and T conjugate channels are 
obtained from d23)-(B0|) by the exchange 8 — ► —8. 



III. 



DERIVATION OF EXACT FORMULA FOR 



In this section, we present three identities on the prod- 
uct of the MNS matrix elements. Then, we derive an ex- 
act formula for the oscillation probability. Furthermore, 



we confirm that the Naumov-Harrison-Scott identity can 
be also derived from our identities. 

The Hamiltonian in matter H is given by 



H = H 




(31) 



where H is the Hamiltonian in vacuum and is diagonal- 
ized as (|), a= 2^/2G F N e E is the matter potential, Gf 
is the Fermi constant and N e is the electron density in 
matter. Note that the second term in (31) is energy inde- 
pendent because the matter potential a is proportional to 
the energy E. The Hamiltonian in matter is diagonalized 
by the effective MNS matrix as 



1 



Ai 



mm = - 1 \-2 



A, 



(32) 



where Ai is mass eigenvalue in matter. The amplitude 
and probability for v a — > vq transition in matter are 
obtained by the replacements mf — > Ai and U — *■ U in 

As one of the approaches to derive the probability, one 
can calculate single U a i by directly diagonalizing H and 
calculate J^g, which is the product of four Us. However, 
the expressions for the probabilities obtained in this ap- 
proach are complicated and it is not easy to extract the 
information on matter effects and 8. Here, we derive the 
probability using another approach. Note that UaiU^ 
appears in the amplitude in (^) and the probability in 
(|TT|). So, we only have to calculate the product UaiUL 

and the expression for single U a i is not necessary. We 
use the relation ( [Ii"| ) in order to calculate U a iU* Bi . We 
introduce p a g and q a g as 



Pap — 2EH a p, 



(33) 



la/3 



{2E) 2 H a B — (2E) 2 (H 1 gH ai — HagH^iiA) 



where (a/37) = ( e M r )j (A tre )j ( re M)- Then, we obtain 
three kinds of identities on U a iU* Qi . First, from unitarity 
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we obtain 



U a \Uf}i + U a2 U^ 2 + U a3 Up 3 — S a p. 



(35) 



Second, 2EH a p = p a p in (p3|), which is also seen in Ref. 
p9| , is rewritten as 

MUalUp-y + \-lUaflUli + \ 3 U a3 Up 3 = pa/3. (36) 

Third, (2E) 2/ H a p = q aj3 in (H), is rewritten as 

^MUalUpx + + ^l^2Ua?,Up 3 — q a/ 3, (37) 

where we use the relation 

(2E) 2 H = (2E) 2 H- 1 (detH) 

= Udmg LL, i j-} W x A!A 2 A 3 . (38) 



Eqs. (p5|)-(|3"7|) on UdU^ can be simply solved as 

PapK + <ia[l — S a p\i(\j + Afc) 



UaiUpi — 



A™ Afc^ 



(39) 



where (ijk) takes (123), (231), (312). Note that p a p and 
Qa/3 become constants in the following cases, 

p a (3 = Pap (a e or f3 e), (40) 
q a /3 = qa/3 (a = e or /3 = e) (41) 

from (|lj). In other cases, H ee is included in p a p and 
q a p, and they depend on a explicitly. Similar expression 
for UaiUpt in ( |39| ) has been obtained by Ohlsson and 
Sncllman in Ref. jl6| and also by Harrison and Scott 
||30f although the method of derivation is different. We 
have given an exact and simple formula of P(y e — > v^) 
for the first time in Ref. p3{ . Here, we give the general 
formula of oscillation probability P(y a — > z/g) applicable 

to all channels by calculating ReJ^g and J from J^g = 

U a iUpi(UajUpj)* . In the case of a 7^ /3, the expression is 
given by 



P{Va 



where 



lie r 3 



cyclic 

E 

m 



v fi) = - 4 E ReJ^sin 2 A^.±2 ]T Jsin2A^ 



cyclic 

E 

m 



(42) 



|Pa/3| 2 AiAj + |g Q/3 | 2 + Re(p a/3 q ; ; / 3)(A l + Aj) 



A l? A 12 A 23 A 3 i 



j = Im(p e pge M ) 
A12 A23 A31 



A' - A ^' £ 
« - 4 £ ' 



(43) 
(44) 



The concrete expression for Ai [B5| is calculated by the 
MNS matrix elements and the masses in vacuum as 



Ai 
A 2 
A3 



-\Vs 2 -3t u+y/3{l-u 2 ) 



y/I 2 



3t 



x/3(l - u 2 ) 



3i, 



(45) 
(46) 
(47) 



where s,t,u are given by 

s = A 2 i+A 3 i+a, (48) 
* = A 21 A 3 i + a[A 21 (l - s 2 2 c 2 3 ) + A 31 (l - s 2 3 )], (49) 



II = cos 



1 _ x f 2s 3 - 9st + 27aA 21 A 31 c 2 2 c 2 3 
3 C0S I, 2(s 2 - 3i) 3 / 2 



(50) 



in the standard parametrization |T^, [li], [[5] . One can see 
that the effective masses in matter do not depend on the 
CP phase 6. 

Finally, we reproduce the Naumov-Harrison-Scott 
identity from (El) as 



A 12 A 23 A 31 J = 



1 



-lxa[H e[l (H eT H T n — H etl H TT ) 



(2Ef 
1 

(2£) 3 ""Peptf^rel 
1 



A 12 A 23 A 3 iJ. 



(51) 



Thus, the Naumov-Harrison-Scott identity is easily de- 
rived from the identities which we used. 



IV. CP DEPENDENCES OF OSCILLATION 
PROBABILITIES IN MATTER 

In this section, we calculate the oscillation probabilities 
in all channels by using (^)-(^J) and investigate their CP 
dependences. 



Let us consider the CP dependence of P{v e — ► v^) 
which is the best channel for the observation of CP vi- 
olation. First, we calculate p efl and q &>1 from ( |36"| ) and 
( p7| ) . As p eiJl and q efl do not include H ee , they are equiv- 
alent to the constants p efl and q e(1 , which are represented 
by the quantities in vacuum. Therefore, we obtain 

Pe)i = A 21 c/ e2 c/; 2 + A 31 u e3 u; 3 , (52) 

q e „ = A 3 iA 2 iC/ e iC/; i (53) 

from (|36l ) and (|37|). Next, we divide p e/J and q eil into the 
terms including and not including 5 as 

Pe»=Pe»e +Pe^ 1e» = q ei fi + 9e M ' ( 54 ) 

where p^, p b e)l , q% and q b eu are real numbers given by 



Pe M = ( A 3i - A 2 is 2 2 )s 23 si 3 ci 3 , 

Pe M = A 2 iSi 2 Ci 2 C 23 Ci 3 , 

q% = -A 3 iA 2 ic 2 2 s 23 si 3 ci 3 , 

9e M = -A 3 iA 2 iSi 2 Ci 2 C 23 Ci 3 , 



(55) 
(56) 
(57) 
(58) 



5 



in the standard parametrization. 

The oscillation probability can be written in the form 



P(v e — > v^) = A eii cos S + B sin 5 + C, 



(59) 



from (fl2|)-(j44|). Note that P(v e — > v^) is exactly linear 
in sin S and cosi5. 

Next, we give the exact expression for A efl , B and C eM . 
We rewrite them in the product of the oscillation part 
which depends on L and (A r )ij, B r and (C r )ij as 



cyclic 

ie M = 51 (ir)ySin 2 A^, 

(tO 
cyclic 

B = J2 Brsm2A^ 

cyclic 

(to 



(60) 
(61) 
(62) 



B is expressed in the form of the sum as (61). Under the 
condition x + y + z = 0, the relation 

sin 2x + sin 2y + sin 2z = — 4 sin x sin ?/ sin z (63) 



holds and -B from ( |si"| ) is rewritten in the form of product 



as 



cyclic 



B = J B r sin2A- J = sin A' 12 sin A' 23 sin A 31 (64) 



A etl is rewritten in the same way. Under the same con- 
dition as in deriving B, the relation 

sin 2 x = — (sin a; sin y cos z + sin x cos y sin z) (65) 



holds and A efi is rewritten as 



cyclic 

A efl = (^)tf sin2 Ay- 
W) 

cyclic 

= - ^r)jk + (A r )ki] cos A^ sin A^ sin A' fci .(66) 

(yfc) 

Substituting (§§)-(§§ into p eM and q ef , in @-(@), I ejU , 
_B and C e ^ are rewritten with the masses and mixings as 



Ae^j, — 



B = 



cyclic 



5[J r A2iA 3 iA fc (A fc - A 3 i) + (i eM )fc] 



(ijk) 

8J r A 12 A 23 A 31 
Ai2A 23 A 3 i 



cosA-, sinA^ sinA' fei , 



sin A' 12 sin A 23 sin A 31 , 



cyclic 

£- 

(y) 



4[A 31 S 13 S 23 C 13 AiAj + (Cep)ij] -2^' 



A r ,A 12 A 23 A 31 



(67) 
(68) 
(69) 



See Appendix for the expression of (A e ^)k and {C efi )ij. 
Note that these expressions (p7j)-(|69"|) are still exact. In 
the limit of small A 2 i, (A efJ- )k and (C e ^)ij are higher 
order in A 2 i and can be ignored. 

Finally, we obtain the well known approximate formula 
by neglecting the smallest effective mass. Considering the 
energy of neutrino (E > 10 GeV) in neutrino factory ex- 
periment and the earth matter density (p ~ 2.8g/cm 3 ), 
matter potential is given by 



= 2V2G F N e E = 7.56 x 10~ 5 eV 2 



P 



E 



gem 3 GeV 



> 2.1 x 10 _3 eV 2 . 



(70) 



It means that the smallest effective mass Ai is almost 
equivalent to A 21 . Other effective masses A 2 and A 3 , 
correspond to a or A 3 i. Accordingly, the coefficients are 



approximated by 



A F 



8J r A 21 A 31 
a(A 3 i - a) 



cos A 31 sina'sin(A 31 — a)', (71) 



B 



8J r A 2 iA 3 i . , . , . , 

—- -smA 31 sma sm(A 3 i -a) , (72) 

a(A 3 i - a) 



a 



Cfi 



A A2 2 2 2 

—— sin (A 3 i - a) . 

(A 3 i - ay 



(73) 



Although the approximate formula derived here is in 
agreement with the ones seen in Ref. J?], ^] , the derivation 
is rather simple. The coefficients A efi and B responsible 
for S are <3(A 2 iSi3). On the other hand, C efl is 0(s 2 3 ). 
Therefore, the terms related to S become relatively large 
in the case of small S13. 
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P(v e -► v T ) 

In this subsection, we calculate P(i> e — * f T ) in the 
same way as P(y e — > ^ M ). We need to learn UeiU^ for 
this purpose. 

At first, we consider p eT and q eT . Since p eT and <7 eT do 
not include H eei they are equivalent to the constants p eT 
and q eT . Therefore, we obtain 

PeT = A 21 u e2 u; 2 + A 31 u e3 u; 3 , (74) 

q eT = A 31 A 21 U el U* Tl (75) 



As in the case of P{v e 
mate formula as 



Vu), we obtain the approxi- 



4A 



31 L 23 J 13°13 
(A 3 1 - ^ 2 



sin 2 (A 31 - a)', 



(83) 



by neglecting the smallest effective mass Ai ~ A21. As 
the result, the leading term C eT which is not related to 
5 is 0(s 2 3 ) as in the case of P{y e — ► v^). It means that 
CP effect becomes relatively large also in this channel. 



from (|36| ) and (|37|). Next, we divide p er and q er into the 
terms including and not including 5 as 



P, 



Per, 



<2e 



(76) 



where Pg T , Pg T , g",. and q\ T are real numbers given by 



Per = ( A 31 - A 2 lS 12 )c 2 3Sl3Cl3, 
Per = -A21S12C12S23C13, 

q~ T - ■ ■ 1 



-A 31 A 2 icl 2 c 23 s 13 c 13 , 
A31A21S12C12S23C13, 



(77) 
(78) 
(79) 
(80) 



in the standard parametrization. Note that p eT and q eT 
are obtained by the replacements S23 — > c 23 , c 23 — > — S23 
in p eM and q efl . The coefficients of cos 5 and sin 5 have the 
opposite sign by the replacements, from (|67j ) and (68). 
Therefore, P(f e — > f T ) can be also written in the form 



In this subsection, we calculate P(vn ^ v t)- The ex- 
pression for the probability of this channel is useful to 
analyze the appearance experiments of MINOS |H| and 
CNGS [[32). It is also needed for the analysis of neu- 
trino factory experiments in which v T is produced by the 
high energy beam. In order to derive Piy^ — ► u T ), 
we calculate U^ll*^ Let us first consider p^ T and q^ T in 
UftiU*^ Since q^ T includes H ee , it is not constant and 
depends on the matter potential a explicitly as 



f 



(225) 



5e) h » t ) 



(2EJ 



;{qp,r - ap^r)- 



(84) 



P(y e -> v T ) = -A efl cos (5- S sin 5 + C eT , (81) Then, is given by 



where 

cyclic 
C er = ^ 

«) 



— 4[A3 1 S 2 3 C2 3 C 2 3 AiAj + (Cer)tjJ \' 



A i7 A 12 A2 3 A 3 i 



sin 2 A' tf (82) 



and the concrete expression for (C eT )jj is given in Ap- 
pendix. As for A e/J and £?, see ( |67j ) and (picf). Note that 
this expression is also exact. In the limit of small A21, 
{Cer)ij is high order in A21 and can be neglected. 



U 77*. = 



PprAi + q^ T _ Ppr(Ai - a) + q^ T 



(85) 



Here, the different point between P{v e — > v^) and 
P{v e — > v T ) is that the coefficient of p^ T is not Ai but 
Ai — a. Fortunately, the matter potential a does not de- 
pend on 5, so we have only to investigate p^ T and q^ T , in 
order to learn the CP dependence of the probability. In 
the standard parametrization, we obtain 



,-«5 



Put = P^e ~ +p" flT +p c fiT e iS , q^ = g£ T e 10 + g£ T + g£ T e'°, 



,-«5 



(86) 



where 



vlr = -A 2 iSi 2 ci2C2 3 si 3 , p\ T = [A 3i c 2 3 - A 2 i(c 2 2 - s 2 2 s 2 3 )]s 23 C2 3 , p c ^ T = A 21 si 2 a 2 sl 3 s 13 , (87) 
s£t = A 3 iA 2 iSi 2 ci2C2 3 si 3 , q^ T = A 31 A 21 (-s 2 12 + c 2 12 s 2 13 )s 23 c 23 , q° T = -A 3i A 2 iSi 2 ci2S2 3 si 3 . (88) 



And then, the transition probability can be written in the form 



Pfir = A^ T cos 5 + B sin S + C^ T + D cos 25, 



(89) 
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where we obtain the expression for the coefficients after some straightforward calculations as 

V = if - 8 ^ A2lA ^ - ^ - a - ^3 - gjs) + OVk] C o S A^ . sin A' sin A^, (90) 
^ — ' A A 

(ijfc) i fe fci 

£ = y^ C -4[A§ 1 a| 3 ^ 3 cf 3 (Ai - o)(Aj - a) + ((Vjjjj 2 ^, 
MT ^ A A 12 A 23 A 31 

n -8A|i(Afc - a - A 3 i) 2 s 2 2 c 2 2 g 2 3c 2 3 s 2 3 -, ~, . -, 

£> = 2^ cosA y -smA it smA K , (92) 

(ijk) A jfc A fei 



and the concrete expressions for (A^ T )k and {C^ T )ij are 
given in Appendix. As for B, see (|68|). Note that these 
expressions (90)-(^) are also exact. In the limit of small 
A 2 i, (A MT )fe and (C eT )ij are high order in A 2 i and can 
be neglected. 

Here, let us comment on the magnitude of each term. 
D is 0(A 21 Si3) and hence it is difficult to observe the 
term proportional to cos 26. A^ T is suppressed by A 2 iSi3 
as in the previous cases. In addition, is also propor- 
tional to c 23 — s 23 . Since the 2-3 mixing is almost max- 
imal from the atmospheric neutrino experiments, A^ lT is 
considered to be rather small. On the other hand, 
has not any suppression factor. Accordingly, the lead- 
ing contribution comes from C^r and the next to leading 
comes from B. Other contributions from A^ T and D are 
extremely small and we give the approximate formula for 

a 



flT- 



As in the previous cases, we obtain 



Cut — 4s 23 c 23 c 13 sin Ao 



i 3i> 



(93) 



by neglecting the smallest effective mass Ai ~ A 2 i. We 
conclude that it is considerably difficult to observe the 
CP effects in the standard scenario. 



i> e ), P(f M — ► v^) and P(i>t — ► v-r) 



In this subsection, we calculate survival probabilities 
using unitarity and the transition probabilities derived in 
the previous subsections. 

Let us first consider P(v e — ► v e ). The expression for 
the probability of this channel is needed, for example in 
neutrino factory experiments with the high energy f e (P e ) 
beam. In the standard scenario, the unitarity relation 



disappear in (p4[). Therefore, P{v e — * v e ) is completely 
independent of 6 and can be written in the form 



P(v e — > v e ) = C ee , 
as first pointed out in Ref. p8fl. Here, 



Cee — 1 C, 
cyclic 

= i+E 



(95) 



(96) 



4[A 31 S 13 C 13 AiAj + (Cee)ij] j , , 



A iJ -Ai 2 A23A 3 i 



sin" A^.(97) 



As for the concrete expression of (C ee )ij, see Appendix. 
Note that this expression is exact. In the limit of small 
A 2 i, (Cee)ij are higher order in A 2 i and can be ignored. 
Furthermore, we obtain an approximate formula 



i 



4 A2 2 2 

(A 3 i - a) 2 



sin 2 (A 31 - a)', 



(98) 



by neglecting the smallest effective mass Ai ~ A 2 i. 
P(y e — > v e ) mainly depends on S13. This channel may 
be interesting because S13 is determined regardless of the 
ambiguity of 6. 



Next, we calculate P(f p 



relation P{v^ 



From unitarity, the 
1 - P(i>n -> v e ) - P(y^ -> v T ) 



holds. Then, the term proportional to sin<5 disappears 
and the probability can be written in the form 

P(v„ -> v») = A^ cos 5 + cos 26, (99) 

where the coefficients are given by 



A„ 



A —A 



-A„ 



(100) 



Cfifi — 1 C e ^ — C^ T — 1 C^ T . (101) 



In the approximation (100) we use the fact that A^ T is 



P(v e ->v e ) = l- P(u e -> Vfj) - P(y e -> v T ) (94) proportional to c^ 3 - sfeirom fl9(J) and is very small in 

t he c ase of maximal 2-3 mixing. In the approximation 



holds. From (|67|) and (|68|) the coefficients of cos 6 and 
sine? are proportional to J r which includes s 2 3C 2 3. As 
P(y e — > v T ) is obtained from P{y e — > v^) by the ex- 



(101), we use the fact that C efl is suppressed by s 2 3 . As 



change s 23 -> c 23 ,c 23 -> -s 23 (i.e. J r 



understand that the terms proportional to cos 6 and sin <5 



for D, it is equivalent to that in Piy^ — > v T ). 

Finally, it may be useful to study the CP dependence 
of P{v T — > v T ) although this channel is very difficult to 
observe. From unitarity, the relation P[v T — > v T ) = 



8 



1 — P{v e — > v T ) — P{Vfj, — ► f T ) also holds. Then, the 
probability can be written in the form 



P(y T — > v T ) = A TT cos <5 + C TT — D cos 28. 
The coefficients are given by 



~A eT - 

1-C e 



A^ T 



-A e 



Cf^T — 1 CflT ■ 



Here, Z) is also the same as that in P(V M — > u T ). 



(102) 



(103) 
(104) 



A uu cos 5 + C u 



P(y T — ► i/ T ) = A TT cos 8 + C T 



- D cos 28, 
D cos 26. 



(109) 
(110) 



Comparing these results with (|25|)-(|30|), we find that the 
probability for each channel in matter has the same form 
with respect to the CP phase 8 as in vacuum. It means 
that the matter effects just modify the coefficients. Fi- 
nally, we comment that the p roba bi litie s for T conjugate 
channels are obtained from (10E)-(11C) by the replace- 
ment 8 — > —8 and those for CP conjugate channels are 
obtained by 6 — > —8 and a — > —a. 



Summary of CP dependences in Matter 

Here, we summarize the CP dependences of the oscil- 
lation probabilities for all channels in matter. They are 
given by 

P{v e -> Vp) = A efl cos 8 + Bsin8 + C efi , (105) 
P(v e -> v T ) = -A^ cos 5 - Bs'mS + C eT , (106) 
P(^ p v T ) = A^ cos 8 + B sin 8 + C„ T + D cos 25(107) 



EFFECTIVE MIXING ANGLES AND 
EFFECTIVE CP PHASE 



In this section, we show that our results ( |105| )-(11C) 
with respect to the CP dependences are reproduced from 
the effective mixing angles and the effective CP phase 
calculated by Zaglauer and Schwarzer |Q. They have 
derived the effective mixing angles and CP phase from 



(108) \U a p\ 



as 



sin 0i3 
sin 2 9 12 
sin 2 023 



A§ - aA 3 + (3 



A13A23 



-(Al 



a\ 2 +/3)A 3 i 



(A 2 - aAi + 0) A 32 - (A 2 - aA 2 + 0)A 3 i ' 
G 2 s| 3 + F 2 c 2 23 + 2GFs 23 c 23 cos 8 



2„2 
*23 



G 2 + F 2 

{G 2 e~ lS 



F 2 e^)s 23 c 23 + GF(c 







F 2 c 2 3 + 2GFs 23 c 23 cos 8) (G 2 c 2 3 + F 2 s 2 3 - 2GF S23 c 23 cos 8) ' 



(111) 

(112) 

(113) 
(114) 



where a, (3, G and F are given by 



a 

P 

Gs 23 
Fc 23 



P\^ + tf 



6 

e ji ' 



(115) 
(116) 
(117) 
(118) 



Let us here describe the derivation of the effective mixing 
angles |i6| . We can easily calculate the 1-3 mixing from 
\U e3 \ 2 = s 2 3 . We can also calculate the 1-2 mixing by 
eliminating £13 from \U e i\ 2 = c 2 2 c 2 3 and |t/ e2 | 2 = s 2 2 c 2 3 . 
Furthermore, the 2-3 mixing can be calculated by elimi- 
nating C13 from the relations, 



\Ur 3 \ 2 



\u e3 u* 3 \ 



\Pe 



\U, 



e?,\ 



A 3 iA 32 (A 2 -aA 3 + /3) 



c 23 c 13- 



(120) 

Next, we consider the CP dependence. One may think 
that the CP dependence of (114) is very complicated due 
to cos 8 in the denominator at a glance. However, our 
results (105)-(110) are very simple and the probabilities 
are linear in sin 8, cos 8 and cos 28. It suggests that the 
complicated CP dependence in 



,-t<5 



cancels with that 
in other terms. From (111)-(114), 8 is only included in 

s'm9 23 except for e~ lS . Furthermore, the numerator of 
sin 2 3 is partially in agreement with the denominator of 



The product of s 2 3, c 2 3 and e lS is calculated as 



M3| 



l^e3^ 3 | S 



be/xA 3 + <? e/ x| 2 



\u< 



c:3 



A3iA 32 (A|-aA 3 +/3) 



~1 ~2 
s 23 c 13' 



(119) 



S 2 3C 2 3e 



-iS 



(G 2 



-iS 



F 2 e td )s 23 c 23 + GF(c 



'23 ^ 



G 2 + F 2 



(121) 
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Note that the CP phase in the denominator of e lS com- 



pletely disappears. The imaginary part of (121) is known 
as the Toshev identity [E7f 



tity ( pl| ) and second is the Toshev identity (122). We 
obtain a new identity 



S23C23 sin 5 = S23C23 sin 5, 



(122) 



which is independent of matter effects. 

Finally, we demonstrate how the CP dependences of 
the oscillation probabilities become simple taking the 
channel P{v e — > v^) as an example. In general, the os- 
cillation probabilities in matter are also obtained from 
those in vacuum by the exch ang e 0y — * 6ij and 5 — ► 6. 
Namely, from ( 15 )-(p~7|) and (pi|), P{v e — > v^) in matter 
is constituted by 



12 



ReJ, 



^j 2 4 = 

ReJ% = 
J = 



— (cf 2 — s\ 2 )Jr COS (5 

1-2 ~2 ~2 l~2 ~2 ~2 \ 
+ s 12 c 12 c 13l s 23 s 13 — C 23J 

J r cos S — s~.,» 

— J r cos 5 — c 
J r sin S 



' -23 
> 

-3) 

^2 ~2 -2 -2 
12 i 23 A 13 c 13! 



S 2 s2 ~2 ~2 
12 6 23*13 c 13! 



(123) 
(124) 
(125) 
(126) 



We can see that sin S (cos 6) always appears together with 
J r in Re J*^ or in J. On the other hand, the CP depen- 

de nce fo r the product S23C23e _l<5 become simple as shown 
in (121). Therefore, the CP dependence of the probabil- 
ity is simplified although that of the effective CP phase 
is complicated. Then, we obtain the same results 



A efi cos S 



B sin S + C, 



(127) 



as (p9|). As for other channels, we obtain the same results 
as those calculated in our approach. This is guaranteed 
by the Toshev identity. 



VI. NEW MATTER INVARIANT IDENTITY 



A12A23A31 Sl 2 Ci 2 Sl3C?3 = Ai 2 A23A 3 iSi2Ci2Sl3C?3 

(128) 

related to the 1-2 mixing and the 1-3 mixing by divid- 
ing the Naumov-Harrison-Scott identity by the Toshev 
identity. 



Here, we give another proof of this new i dent ity (128). 
We can independently derive the identity (128) without 
using other identities in this proof. At first, taking ki, fc 2 
as the eigenvalues of sub-matrix 



1 ( H^t 
2E V Htv. H tt 



we rewrite \U e i\ 2 obtained from (B9) as 



\U,,.\ 2 



- (?W + Prr)\ + ge 

AoiAfc, 
(Aj - fciKAj - k 2 ) 



Then, we obtain 



(129) 



(130) 



(A 12 A 23 A 3 l) 2 |C/ e l| 2 |C/e 2 | 2 |C/e3| 2 

= n ( A i- fc o(A 2 -^)(A 3 -fc l ) 

i=l,2 

= Y\_ ( m i — ^«)( m 2 — — ki) = const., (131) 

i=l,2 



In the previous section, we rederived two matter invari- 
ant identities. One is the Naumov-Harisson-Scott iden- 



where we use the relations obtained from the character- 
istic equation of Hamiltonian 



Ai + A2 + A3 
A1A2 + A2A3 + A3A1 
A1A2A3 



„„2 2 , „,2 2 , „„2 2 

m 1 m 2 + m 2 m 3 + m 3 m 1 
rn\rn\rn\ + ak\k2, 



a(ki + k 2 ), 



(132) 
(133) 
(134) 



from the first line to the second line in ( 131 ). Note that 
ki is the eigenvalue of the sub-matrix of Hamiltonian 
in vacuum and then does not depend on matter effects. 
Therefore, (131) is a matter invariant quantity and it is 



rewritten as ( 12S ) 



VII. SUMMARY 



We summarize the results obtained in this paper. We 
consider the probabilities and the CP dependences for all 
channels within the context of standard three neutrino 
oscillations in constant matter. 
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(i) We have derived an exact and simple formula for 
the oscillation probability applicable to all chan- 
nels. We have also investigated the CP depen- 
dences of various probabilities using this formula. 
As the results, P(y e — + v^) can be written in the 
form A efl cos 6 + B sin (5 + C eu . Piy e — * Vt) has the 
same form as P(y e — > v^). On the other hand, the 
probability for to v T transition can be written as 
P(v u — > v r ) = A^ T cos S + B sin S + C^ T + D cos 25. 
In survival probabilities P(y e — > v e ), P{v^ — » iO 
and P(y r — ► i/ T ), there is no sin<5 term. As for 
P(z/ e — > z/ e ), cos 6 term also disappears and the 
probability is completely independent of 5. We 
have found that the probability for each channel in 
matter has the same form with respect to the CP 
phase as in vacuum. That is, the matter effects just 
modify the coefficients A, B, C and D. We have 
also given the exact expression for the coefficients 
in constant matter. 



(ii) We have shown that our results with respect to the 
CP dependences are reproduced from the effective 
mixing angles and CP phase derived by Zaglauer 
and Schwarzer. Although the effective CP phase 
are complicated, the CP dependence of probabil- 
ity is simplified for all channels. This is guaran- 
teed by Toshev identity. Finally we have obtained 
a new identity related to the 1-2 mixing and the 
1-3 mixing by dividing the Naumov-Harrison-Scott 
identity by the Toshev identity. 
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APPENDIX A: SUBLEASING TERMS IN THE 

Here, we describe the terms which include higher order 
of A21 explicitly in A a p and C Q ^. 
First, the terms in P(y e — > v^) are 



{A efl ) k = A 2 21 J r x [A 31 A fc (c 2 2 - sj 2 ) + X 2 k s 2 12 - Aj lC j 2 ] , (Al) 

(Cepkj = A21S13 x [A 31 {~X l (X j sj 2 + A 3l c 2 2 ) - Aj(Ajsf 2 + A 3 ic 2 2 )}s 23 c 2 3 ] 

+ A 2 21 x [(Aj - A 3 i)(A j - A 31 ) s 2 12 c 2 2 c 2 3 c 2 3 ] 

+ A| lSl3 x [(X lS 2 2 + A 31 c 2 2 ){X 3 s 2 2 + A 3l c 2 2 )s 2 23 c 2 3 ] , (A2) 



where the second line of {C e ^)ij is relatively large in the case of A21/A31 < s\ 3 so that there is no suppression factor 

*13- 

Second, the terms related to P(v e — ► v T ) are 

(Cer)ij = A 21 sf 3 x [A 31 {-X l (X j sj 2 + A 31 cj 2 ) - Xj(Xisl 2 + A 31 cl 2 )}cj 3 cl 3 ] 
+ A 2 21 x [(A, - A 3 i)(A J - A 31 ) s 2 2 c 2 2S 2 3 c 2 3 ] 

+ A 2 lS 2 3 x [(AiS 2 2 + A 31 c 2 2 )(A J -s 2 2 + A 31 c 2 2 )c 2 3 c 2 3 ] , (A3) 

where there is also no suppression factor s 2 3 in the second line. 
Third, the terms related to Piy^ — > v T ) are 

{Anr)k = -A2!S 13 (c| 3 - sj 3 ) 

x [si2Ci 2 s 23 c 23 (A fc - a - A 3 i){A 3 i(s 2 2 - s 2 3 c 2 2 ) + (A*. - a - A 3 i)c 2 2 }] (A4) 
(C^ T )ij = A 2 i x [A 3 i{A 3l (Aj + Xj - 2a)(c 2 3 + c 2 2 ) - 2(A, - a)(Xj - a)c 2 2 }s 2 23 c 2 23 c 2 13 ] 
+ A 2 is 2 3 x [A 3 i{— A 3 i(A^ + Xj - 2a) + 2(A, - a)(Xj - a)}s 2 12 s 2 23 c 2 23 cj 3 ] 
+ A 2 ! x [{A 31 c 2 3 + (Aj - a - A 31 )c 2 2 }{A 31 c 2 3 + (A, - a - A 31 )c 2 2 }sl 3 cf. 3 ] 
+ A^sf 3 x [-A 3 i(Aj + Xj - 2a - 2A 31 )s 2 l2 s 2 23 c 2 23 c\ 3 

+ (A, - a - A 31 )(Aj- - a - A 31 ){s\ 2 c 2 l2 {c 2 23 - s 2 23 ) + s\ 2 s 2 23 c 2 23 s 2 13 }\ . (A5) 

Finally, the terms related to P{v e — ► v e ) are 

(C ee )y = A 21 s 2 13 x [A 31 {-Xi(X jS j 2 + A 31 c 2 12 )- X j (X i sj 2 + A 31 c 2 12 )}c 2 13 ] 
+ A 21 x [(A, - A 31 )(Aj - A 31 )s 2 2 c 2 2 c 2 3 ] 

+ A 2 lS 2 3 x [(A lS 2 2 + A 31 c 2 12 ){X 3 s 2 2 + A 31 c 2 2 )c 2 3 ] , (A6) 

I 

where there is no suppression factor s 2 3 in the second line as in P(v e — > v^) or P(y e — > v T ), Therefore, the 
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contribution is relatively large in the case of A21/A31 < S13. 



